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■ Abstract 



We present the calculation of the order corrections to the coefficient functions 
contributing to the longitudinal {Fi{x,Q'^)) and transverse fragmentation functions 
{Ft{x, Q"^)) measured in electron-positron annihilation. The effect of these higher or- 
der QCD corrections on the behaviour of the fragmentation functions and the corre- 
sponding longitudinal {daL^x^Q"^) / dx) and transverse cross sections {daT{x,Q'^)/dx) 
are studied. In particular we investigate the dependence of the above quantities on 
the mass factorization scale (M) and the various parameterizations chosen for the 
parton fragmentation densities D^i^x^M"^) {p = q, g; H = n"^ , K"^ , P, P) . Our anal- 
ysis reveals that the order contributions to Fl{x,Q'^) are large whereas these 
contributions to Ft{x,Q'^) are small. From the above fragmentation functions one 
can also compute the integrated cross sections and in an independent way. The 



sum (7tot = + ctt, corrected up to order a^, agrees with the well known result in 
the literature providing us with an independent check an our calculations. 



1 Introduction 



Semi-leptonic processes represented by electron-positron annihilation into hadrons, 
deep inelastic lepton-hadron scattering and the Drell-Yan process have provided us 
with the most valuable testing grounds for perturbative quantum chromodynamics 
(QCD). Perturbative calculations in next-to-leading order, and in some cases even 
to higher order, give a good explanation of numerous quantities measured in various 
experiments The reason for these successes originates from the experimental as 
well as theoretical characteristics of the above reactions. From the experimental view- 
point semi-leptonic reactions provide us with an overwhelming amount of data and 
in the case of electron-positron annihilation into hadrons and deep inelastic lepton- 
hadron scattering the background is fully under control. Therefore the systematical 
and statistical errors are very small. From the theoretical viewpoint we want to 
mention the following features. First, the Born approximation to semi-leptonic cross 
sections is of purely electroweak origin so that it is independent of the strong cou- 
pling constant a^. Since the electroweak standard model is tested up to about a few 
promille by the LEPl-experiments p| each deviation from the Born approximation 
is due to the strong interactions. Second, if one limits oneself to the computation 
of semi-inclusive or inclusive quantities, like structure functions or total cross sec- 
tions, the final hadronic state is completely integrated over and we do not have to 
care about problems as jet definition or hadronization effects. The third feature is 
that it is possible to extend the calculation of the QCD corrections to the above 
integrated quantities beyond next-to-leading order. Examples are the order con- 
tributions to the coefficient functions corresponding to the Drell-Yan cross section 
da/dQ"^ and the deep inelastic structure functions Fk{x,Q'^) P where denotes 
the virtuality of the electroweak vector bosons 7, Z, W. Order corrections are 
even known for sum rules dx x^~^ Fk{x , Q"^) {n < 10) and the total cross sec- 
tion crtot(e'^e~ — > "hadrons") The reason that these higher order corrections are 
much easier to compute than those encountered in e.g. hadron-hadron collisions (ex- 
cept for the Drell-Yan process) can be attributed to the simplicity of the phase space 
integrals and the virtual corrections appearing in semi-leptonic processes. Moreover 
if one integrates in the latter processes over the total hadronic state one can use 
alternative methods to evaluate the Feynman diagrams (see e.g. |]TB||), which are not 
applicable to hadron-hadron reactions or to more exclusive semi-leptonic processes. 
In the past the order contributions to the coefficient functions have been cal- 
culated for the Drell-Yan cross section da/dQ"^ and the deep inelastic structure 
functions Fk{x,Q'^) 0. However the same corrections were not computed for the 
fragmentation functions showing up in the process e^e~ H + "X" where H is 
the detected hadron {H = tt^,K^,P,P) and "X" stands for any inclusive hadronic 
state. These corrections are needed because of the large amount of data which have 
collected over the past ten years. The above process has been studied over a wide 
range of energies of many different e+e~-colliders. Data have been collected from 
DASP (yi = 5.2 GeV) 1T|, ARGUS (^i = 10 GeV) 0, TASSO (v^ = 22,35,45 
GeV) [^, MARK II ^ and TPC/27 (v^ = 29 GeV) CELLO (v^ = 35 GeV) 



2§, AMY (v^ = 55 GeV) ^ and the LEP experiments DELPHI p^, ALEPH 
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OPAL |3y, |3^ {y/s = 91.2 GeV). In particular the last two experiments found a dis- 



crepancy between the measured longitudinal fragmentation function Fl{x,Q'^) and 
its theoretical prediction computed up to order as- We want to fill in this gap in our 
knowledge by presenting the order a1 contributions to the longitudinal {Fl^x.Q"^)) 
and transverse fragmentation {Ft{x,Q'^)) functions and discuss their phenomeno- 
logical implications. The coefficient functions corrected up to order are already 
computed for Fl{x.,Q'^) and can be found in recent work [^. Here we want to add 
the order a1 contributions to the coefficient functions corresponding to Q^) 
which are much more complicate. The order corrections to the asymmetric frag- 
mentation function will be postponed to a future publication. Although a complete 
next-to-next-to-leading (NNLO) order analysis of the transverse (and also asymmet- 
ric) fragmentation function is not possible, since we do not know the three-loop order 
timelike DGLAP [S^ splitting functions, one can still study the effect of the order 
corrected coefficient functions. Furthermore one can obtain the transverse cross 
section ctt for which analysis the DGLAP splitting functions are not needed so that 
the former is factorization scheme independent. The sum of the transverse (ar) and 
the longitudinal (o"l) cross sections yield crtot(e^e^ "hadrons"). It turns out that 
the (Jtot presented in this paper is in agreement with the order a1 corrected result 
quoted in the literature providing us with a very strong check on our calculations. 
This paper will be organized as follows. In section 2 we introduce our notations of the 
fragmentation functions and the corresponding cross sections. In section 3 we give 
an outline of the calculations of the parton subprocesses contributing to the process 
e^e~ H + "X" up to order a^. In section 4 we perform the renormalization and 
mass factorization of the partonic quantities providing us with the longitudinal and 
transverse coefficient functions. The discussion of our results will be presented in 
section 5 and a comparison with data coming from recent and past experiments on 
electron-positron annihilation will be made. The long expressions obtained for the 
order corrected coefficient functions are presented in the MS-scheme and the A- 
(annihilation) scheme in appendix A and appendix B respectively. 
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2 Single particle inclusive cross sections 



In this paper we want to study the QCD corrections to the single particle inclusive 
process 

e+ + ^ 7, Z + "X" , (2.1) 

where "X" denotes any inclusive final hadronic state and H represents either a specific 
charged outgoing hadron or a sum over all charged hadron species. The unpolarized 
differential cross section of the above process is given by [^0|, ^ 



-(1 + cos^ 0)—!- + - sm^ 9—^ + - cos 9—^. (2.2) 



dxd cos 9 8 dx 4: dx A dx 

The Bj0rken scaling variable x is defined by 
2pq 

Q 



^ = ^' g' = Q'>o, 0<x<l, (2.3) 



where p and q are the four-momenta of the produced particle H and the virtual vector 
boson (7, Z) respectively. In the centre of mass (CM) frame of the electron-positron 
pair the variable x can be interpreted as a fraction of the total CM energy carried 
away by the hadron H. The variable 9 denotes the angle of emission of particle 
H with respect to the electron beam direction in the CM frame. The transverse, 
longitudinal and asymmetric cross sections in ([^.2D are defined by , , and 



respectively. The latter only shows up if the intermediate vector boson is given by 
the Z-boson and is absent in purely electromagnetic annihilation. 
In the QCD improved parton model which describes the production of the parton 
denoted by p and its subsequent fragmentation into hadron H, the cross sections 
[k = T, L, A) can be expressed as follows 



da^ dz 
dx 



p=i ^ 



(2.4) 



for k = T, L. In the case of the asymmetric cross section we have 

^^dz^""' 

z 



dal 
dx 



X 



p=i 



(2.5) 



In the formulae ( p. 41) and ( p.5|) we have introduced the following notations. The func- 
tion D^{z^ M^) denotes the gluon fragmentation density corresponding to the hadron 
of species H. The same notation holds for the quark and anti-quark fragmentation 
densities which are given by Dp{z, M^) and Dp{z, M^) respectively. Further we have 
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defined the singlet (S) and non-singlet (NS, A) combinations of quark fragmentation 
densities. They are given by 



p=i 



(2.6) 



(2.7) 
(2.8) 



The index p stands for the quark species and Uf denotes the number of light flavours. 
Assuming that the charm and the bottom quark can be treated as massless we can 
put Hf = 5 and the indices p = 1,2,3,4,5 stand for p = u,d,s,c,b. Further the 
variable M appearing in D^{z,M^) stands for the mass factorization scale which 
for convenience has been put equal to the renormalization scale. The pointlike cross 
section of the process 



e"^ + e" 



p + p, 



which shows up in (|2.4| ) is equal to 



2 2 I "^Q (Q 



Ml 



\Z{Q 



2^2 



eeCpCv/Cv^p + 



{Q 



2\2 



\Z{Q 



(2.9) 



(2.10) 



p=i 



with 



Z{Q^) = Q^-Ml + tMzTz. 



(2.11] 



(2.12) 



Here A^ stands for the number of colours (A^ = 3) and Mz, T z denote the mass and 
width of the Z-boson respectively. For the latter we have used the narrow width 
approximation. Furthermore we have neglected all quark masses in ( p.lO|) . The 
charges of the lepton and the up and down quarks are given by 



-1, 



(2.13) 



_ 2 

The vector- and axial-vector coupling constants of the Z-boson to the lepton and 
quarks are equal to 



Ca,u 



2 sin26'vK 
—CA,d = —C/ 



-CAiii-^^i^^ew) 



(2.14) 



Cv,u = Cap (1 - 7t sin^ 6w] 



V4 
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where 6w denotes the Weinberg angle. 

The electroweak couphng constants also appear in the asymmetry factor A^'^ ( pl^ 
which is given by 



2g2(Q2-M|) 



+ 4 T) C A pCa vCv eCv 



^ (2.15) 



The QCD corrections in (|2.4|), (|2T5| ) are described by the coefficient functions C^^ 
{k = T,L,A; i = q, g) which can be distinguished with respect to the flavour group 
SUijif) in a singlet (r = S) and a non-singlet part (r = NS). They depend on 
the factorization scale M and in order a1 on the number of flavours nj. As will 
be shown later on the gluonic coefficient function only receives contributions from 
flavour singlet channel partonic subprocesses so that we can drop the superscript S 
on Cg. However the quark coefficient functions can be of flavour singlet as well as 
flavour non-singlet origin. Up to first order in the strong coupling constant as it 
turns out that C^^ = However in higher order both quantities start to deviate 
from each other. Hence we define the purely singlet coefficient function C^^ via 

^l, = ^Z + ^Z- (2-16) 

Like £-k,g the purely singlet coefficient function only receives contributions from the 
flavour singlet channel partonic subprocesses which for the flrst time show up in order 

Using charge conjugation invariance of the strong interactions one can show that 
C^s^ = -C^ll and C^^^ = C^,,, = 0. This implies that to (gj) only non-singlet 



channel partonic subprocesses can contribute. Another important property of the 
coefficient function is that they do not depend on the probe 7 or Z or on the elec- 
troweak couplings given in (|2.13|) , ( p.l4|) so that one can extract the overall pointlike 



cross section crp"-* (|2.10| ) or the asymmetry factor A^'^ ( p.l5|) . However this is only 



true if all quark masses are equal to zero and if one sums over all quark members in 
one family provided the latter appear in the inclusive state of the partonic subprocess 
(see section 4). 

From (|2.2| ) we can derive the total hadronic cross section 

(^') = 2l^X (^ = ^'^'^)' (2.18) 



CTtot 

with 

CTfc 



where one has summed over all types of outgoing hadrons H . Hence we obtain the 
result 

a,M^) = Reea^,{Q^), (2.19) 
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where i?ee represents the QCD corrections to the pointhke total cross section a['^l{Q 



At this moment the perturbation series of Re 
Up to order it reads 



is aheady known up to order 



llln|^-44C(3) 



1 23 

+ — K nfCpTf 



41n|^ + 16C(3)-22 



(2.20) 



In section 4 we also want to present the coefficient functions C^/ up to order and 
show that they lead to the same Ree as calculated in the literature (see section 5). 
Finally we also define the transverse, longitudinal and asymmetric fragmentation 
functions F^{x,Q^)f\ 



(J- 



dx 



k = {T,L,A). 



Further the total fragmentation function is given by 
F''ix,Q')=F^{x,Q')+F^{x,Q'). 



(2.21) 



:2.22l 



In the case the virtual photon dominates the annihilation process ( p^.l| ) one observes 
that, apart from the charge squared in (|2.10|) {p = u,d), the above structure func- 
tions are just the timelike photon analogues of the ones measured in deep inelastic 
electron-proton scattering. When the Z-boson contributes we will define in section 4 
for each combination of the electroweak coupling constants in (|2.10| ) a separate struc- 
ture function. However for the discussion of our results in section 5 this distinction 
will not be needed. 



* Notice that we make a distinction in nomenclature between the fragmentation densities and 
the fragmentation functions . 
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3 Fragmentation coefficient functions in e an- 
nihilation up to order 

In this section we will give an outline of the calculation of the order corrections 
to the fragmentation coefficient functions. The procedure is analogous to the one 
presented for the calculation of the Drell-Yan process in 0] and the deep inelastic 
lepton-hadron reaction in 0. The coefficient functions originate from the following 
reaction 

V{q) ^ "p(fco)" + Pi{h)+P2{k2) + ---+Pe{ke), (3.1) 

where V = j, Z and "p" denotes the detected parton which fragments into the hadron 
H. The process ( |3.1|) is inclusive with respect to the partons pi {i = 1,2,. . . ,i) so 
that one has to integrate over all momenta indicated by ki. Notice that the first part 
of reaction ( p.lD i.e. —>■ V is not relevant for the determination of the coefficient 

function. 

Up to order all parton subprocesses represented by ( p.l|) are listed in table From 
the amplitude M^{i) describing process ( p.lD one obtains the parton structure tensor 
(indicated by a hat) 

n— 3 °° f 

Wl^V'ip, = V 5Z / dPS^'^ M;(£) Mf (£)*. (3.2) 

i=i 

Here / dPS'-'''' denotes the fc-body phase space integral defined by 

IPS^^) = jn / (^^^(^J)} (27r)"5(")(g- A:o-n fc.), (3.3) 

5^{k]) = e{k^)6{k% (3.4) 

and and u stand for the Lorentz indices of the vector bosons V and V respectively 
with 1/ = 7, Z and V' = 7, Z. Further we have defined the partonic scaling variable 

= Q^. (3.5) 

and the factor z^~'^ in ( p. 21) originates from the n-dimensional phase space of the 
detected parton p ( ^TTf) . It appears in the definition of the cross sections dak,p/dz 
which are the partonic analogues of the hadronic cross sections in ( |2.2| ). The former 
are proportional to the functions J^k,p defined below. 

To regularize the ultraviolet (U), infrared(IR) and collinear (C) divergences showing 
up in expression (|3.2| ) we have chosen the method of n-dimensional regularization. 
Therefore the phase space integral in ( p.3| ) is generalized to n dimensions so that 
the above divergences show up as pole terms of the type (l/e)"^ with e = n — 4. 
The calculation of the matrix elements M^{k) My{k)* was performed in n dimensions 
using the algebraic manipulation program FORM After having computed the 
traces we have to integrate the matrix elements over all internal loop and final state 
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momenta where the momentum of the detected parton is kept fixed. In this paper 



we take all partons to be massless. The case of massive quarks is discussed in pO 
where their contributions are presented up to order as- 
The parton structure tensor in ( p.2|) can be also written as 



1 °° r 

w^ir''"Hfco,g) = ^E / dPsW(o|jf)(o)b,{M)(p,{MUrHo)|o), (3.6) 



where jji^'' is the electroweak partonic current corresponding to the vector boson V. 
Using Lorentz covariance and CP invariance (|3.6|) can be written as follows 



91 92 



2kQq 



2{koqy 



■ :FaAz,Q^)- (3.7) 



We will call J^k,p{z, Q"^) {p = q, g) the parton fragmentation functions which describe 
the Born reaction plus the higher order QCD corrections represented by the parton 
subprocesses in table ^ The vector and axial-vector couplings of the quark q interact- 
ing with the vector boson V are given by Vq^^ and a^^^ respectively. In the standard 
model they read 

^,(7) _ 2 



,(7) 



<Ju 



1 

■3' 



i-lsin^^ 



■\ + fsin^ 





= 0, 




= 0, 




_ 1 

2' 







As we have already mentioned in section 2 below ( p.l6|) in the case of massless quarks 
the electroweak factors can be completely factorized out of the radiative corrections 
according to ( |3.71 ) so that jFk,p {k = T, L, A) do not depend on them. Therefore 
we can also put them all equal to 1/a/2 without affecting the parton fragmentation 
functions. Hence the latter are obtained via the following projections 



^T,p{^, Q^) 



n-2 
1 



g2 ^ koq ' ' 



koq 



g2 (n-2)(n-3) 



flUl 



(3.9) 
(3.10) 
(3.11) 
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where according to the prescription in we have contracted the Levi-Civita tensors 
^tj.t^af3 n-dimensions. In this paper we will compute the transverse and longitudinal 
fragmentation functions only and leave the calculation of the asymmetric fragmenta- 
tion to a future publication. The computation of the latter involves the prescription 
of the 75-matrix and the Levi-Civita tensor in n-dimensions which is quite intricate. 
We will now discuss the QCD corrections order by order in perturbation theory. In 
zeroth order in (see fig. ^ we obtain the simple parton model results 



(0) 



T,q 



6{l-z), 



0; 



7(1) 



0. 



(3.12) 



The first order corrections denoted by JF^ / {k = T, L;i = g, g) have been calculated in 



the literature ^ . In the case of n-dimensional regularization they are computed 
up to finite terms in the limit e — > and can be found in Since the mass 

factorization has to be carried out up to order one also needs those terms in 
^l^^{z, Q^, e) which are proportional to e. Therefore we have repeated the calculation 



of the graphs in fig. H and 



(1) 



and the results can be presented in the following form 

e/2 

(3.13) 



e a 



(1) 



e/2 



p(0)l I c^l) I 
^qq ^ ^ ^T,q ^ 



e a 



(1) 

T,9 



(3.14) 



(1) 



(1) 



An 
47r 



e/2 r 



e/2 



gq 



-T,g 



(1) 
\9 



(3.15) 



(3.16) 



The pole terms 1/e stand for the collinear divergence in the final state and /i^ and 5*^ 
are artefacts of n-dimensional regularization. The mass parameter /i originates from 
the dimensionality of the gauge coupling constant in n dimensions and should not be 
confused with the renormalization scale R and the mass factorization scale M. The 
spherical factor Se is defined by 



5^ = exp 



-e{'-/E — ln47r) 



(3.17) 



Further cts denotes the bare coupling constant and Pj;p = q,q,g) stand for the 
lowest order contribution to the DGLAP splitting functions |3^. Using our conven- 
tion they are presented in eqs. (2.13)-(2.16) of 0]. Notice that in lowest order there 
is no difference in the expressions for Pj^p found for the deep inelastic structure func- 
tions (spacelike process) and those appearing in the fragmentation functions (timelike 
process). In next-to leading order the DGLAP splitting functions are different for 
spacelike and timelike processes as will be shown later on. 

The coefficients cj!'], presented in the MS-scheme, are already calculated in the lit- 
erature ^ (see also appendix A). Furthermore we also have to compute the 
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coefficients a[,\- (proportional to e), since they are needed for the mass factorization 
which has to be carried out up to order a^. The results are 

a^l\ = CF {-l + ln(l-2) +21nz}, (3.18) 



= Cf I D,iz) - ^D,iz) + Q - 3C(2) ) Doiz) - 1(1 + z) \n\l - z) 

1 + z"^ 1 _|_ ^2 o 

+ 2- ln2ln(l-;z) + 2- z - 3- \nz + -{l-z) ln(l - z) 

+ 3(1 - ^) In z - ^ + + ^(1 + .)C(2) + 5(1 - ^)(9 - ^C(2)) I , (3.19) 



a^l]^ = Cf \ 4^^(ln(l - z) + 2hiz - 2) \ , (3.20) 



ajlj^ = <! (7 - 2 + z)(ln^(l - z) + 41n2ln(l - z) + z - 3C(2)) 



z 



1- z ^ 

-4 (ln(l-;z)+21n;z-3)+4;z y. (3.21) 

The calculation of the order corrections proceeds in the following way. First we 
have the two-loop corrections to the quark-vector boson vertex represented by the 
graphs in fig. |^ which only contribute to ^^q- '^^^ two-loop vertex correction can 
be found in eq. (2.49) of (see also appendix A of 0]). The result agrees with 
the one quoted in [^]. Notice that the first graph in fig. ^ does not contribute for 
V = J because of Furry's theorem. It only plays a role in the case V = Z provided 
one sums over all fiavours in a quark family in order to cancel the anomaly which 
originates from the triangle fermion sub-loop. Since all quarks are massless the final 
result for this graph is zero too even in the case oiV = Z. 

Next we have to compute the one-loop virtual corrections to the radiative process 
in fig. § which contribute to JF^^^^ as well as jF^^j [k = T,L). The corresponding 
graphs are shown in fig. ^ Notice that we have omitted the diagrams with the 
self energy insertions on the external quark and gluon legs. Their contributions 
vanish because of the method of n-dimensional regularization and the on-mass shell 
conditions = kj = 0. Another vanishing contribution happens for the last graph in 
fig. ^ when V" = 7 because of Furry's theorem. In the case of = Z it only contributes 
when the quarks are massive. However here one has to sum over all members of a 
quark family in order to cancel the anomaly originating from the triangle fermion 
loop. 

The amplitude of the parton subprocesses in fig. |^ will be denoted by M(2) (see 
( p.l|) where i = 2). The momenta of the incoming vector boson V and the outgoing 
partons are parameterized like 

q = ys(l,0„_i), 
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A;i = ^— (I,cos6'i,sin6'i,0„_3), k2 = q - h - ki, (3.22) 

where 0„ stands for the n-dimensional null vector. The phase space integral in (|37^ ), 
( p.3|) becomes 

. f dyv"\l-v)"^ \M(2)\\ (3.23) 

(3.24) 

Here we have defined the Lorentz invariants 

s = Q\ si = {ko + kif, S2 = {ko + fc2)^ si2 = {ki + kif, (3.25) 



with s = Si + S2 + Si2- The parameterization of ( |3.23| ) follows from momentum 
conservation and the on-shell condition k^ = kj = 0. Hence we get 

S2S12-S1S /I ^ /T \ 

cos6'i = - -, si = z{l-y)s, Si2 = {l-z)s, S2 = zys. 

{s - Si2){s - S2) 

(3.26) 

The Feynman integrals corresponding to the one-loop graphs which contribute to 
M(2) in ( |3.23| ) contain loop momenta in the numerator. They can be reduced to 
scalar one-loop integrals using an n-dimensional extension of the reduction program 
in |H9[ . The expressions for these scalar integrals which are valid for all n can be found 



in appendix D of [47|. The phase space integrals which emerge from the computation 



of |M(2)|^ are very numerous so that we cannot present them in this paper. They 



are calculated algebraically using the program FORM 

The most difficult and laborious part of the calculation can be attributed to the 
parton subprocesses ( |3.1| ) where one has to integrate over three partons in the final 
state (see also table p. These parton subprocesses are depicted in figs. ^ providing 
us with the amplitude M(3) (see p.l| ) where i = 3). The graphs in fig. ^ determine 

JF^^^ as well as JF^^j {k = T, L) whereas the graphs in fig. ||, which only contain quarks 

" (2) 

and anti-quarks in the final state, contribute to J-'l ^ only. For the computation of 
the three body phase space integrals we choose the following parameterization for the 
momenta of the virtual vector boson V and the outgoing partons (see P5| ) 



q = a/s(1,0„_i), 

■5023 — S23 



[1,1,0. 



n~2 



11 



h = { ^^l^ ^^•^ ) (l,cosx,sinx,0„_3), 



h = -^^5^(1, cos^i, sin 6*1 cos 6*2, sin sin 6^2, 0„_4), 

h = -^^5^(1, - cos6'i, - sin6li cos6'2, - sin sin ^2, 0„_4), (3.27) 
where we have defined the invariants 

^ij — ij^i ~\~ kj) , Sijm = {ki + kj + k^) , s = q . (3.28) 

From momentum conservation and the on-mass shell conditions one can derive 

2S23(S + S23 — S023 — ■5123) . . 

1 _ COS X = TT (3.29) 

(S023 - S23)(Sl23 " S23) 

The three-body phase space integral in (|3.2|), (p.3|) can be expressed as 



/ 



dPS(3) |M(3)|' = ^ ttA-^ s'^' z^^'I\\ - zy+'- 
' ^ ^' 2%4 T{l + e) {AttY ^ ' 

' dy, f dy,y\+^{l - y.f'^yfil - y^f" (1 - 1/2(1 - z))' 
^0 



£-2 



•/ d9i{sm9iY+' c/^2 (sin ^2)" 1^(3)1% (3.30) 
Jo JO 

where the invariants in (|3.28|) depend on z, yi, and y2 in the following way 

yizs z{l-z)yiy2S 

S023 = z 7Z r, S123 = (1 - Z)S, S23 = 7Z r, 

1-1/2(1-2) 1-1/2(1-2;) 

2koq = zs, soi = z{l - yi)s. (3.31) 



Before we can perform the angular integrations the matrix element |M(3)|^ has to be 
decomposed via partial fractioning in terms which have the general form 

T-inin2n3n4 — (q ^"i I'e "l^a l' e (a "l"* 

— K^iiji) \^i232) K^njsJ \^MJ4) 5 

= ■■■,2,l,0,-l,-2,---. (3.32) 

The decomposition can be done in such a way that one invariant e.g. Sj^j^ in the 
product ( ^.32| ) depends on the polar angle 9i whereas an other invariant e.g. Si^j^ 
contains the polar angle 61 as well as the azimuthal angle 62- The remaining invariants 
i.e. Sjgjg and Si^j^ do not depend on the angles. 

Sometimes it happens that the azimuthal angle 62 also appears in Sj^j^. In this 
case one has to rotate the frame in (|3.27| ) so that Sj^jj becomes independent of the 
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azimuthal angle. This is always possible because the phase space integral ( p.4| ) is 
Lorentz invariant. The angular integrals take the form 

T(^'^)=rdO CdO (singQi+^sing,)- 

io io (« + ^cos^i)H^ + 5cos^^i + Csin^icos^2)^' ^ ' 

where a, 6, A, B, and C are functions of the kinematical invariants s, S123, S023, S23 
( p.31|) . These integrals can be found in appendix C of However they have to be 



extended by including terms which are proportional to e where the degree k has to 
be larger than the one appearing in the integrals of |^|. This is necessary because 



these terms contribute due to the appearance of high power singularities (l/s)^ in the 
phase space integral ( p.30| ). The n-dimensional expression for ( p.33| ) becomes very 



cumbersome if 7^ and A"^ ^ B"^ + C"^. Fortunately this situation can be avoided 
when one chooses the frame presented in ( p.27| ) . In this frame the worst case is given 



by 7^ 6^, = B'^ + C'^ or = b'^ , ^ B'^ + C"^ . These type of integrals have to be 
partially done by hand before one can algebraically evaluate expression ( |3.30D using 
the program FORM ||4^. The angular integrals are easy to perform when = 6^ 



and A"^ = B"^ + because they can be expressed into a hypergeometric function 
2-Fi(a, /?; 7; x) |5^. Inserting the latter in ( p.30| ) the remaining integrations are then 



again performed using the algebraic manipulation program FORM. 
Finally we would like to comment on a special type of term appearing in the transverse 
parton fragmentation function TT,q{z-,Q'^ ■,£)■ They only show up in the non-singlet 
part and the order a™ contribution takes the form 

2m— 1 f / \ 

n':^{z,Q^e)=J2il-z)^^-'l^, (3.34) 
e=-i 

where /^(l) is finite. These type of terms originate from gluon bremsstrahlung 
(figs. §,H,|^) and gluon splitting into a quark-anti-quark pair (fig. |). In the limit 
2; — >■ 1 all gluons become soft and the angle between the quark and anti-quark pair 
goes to zero (collinear emission). In the next section expression (|3.34|) has to be 
convoluted with the so-called bare fragmentation densities D^{z) (for the definition 
see section 4) which yields the integral 



2m— 1 



dzb^r-) (l-z)f-i^. (3.35) 



X 



Inspection of the above integral reveals that aX z = 1 one gets an additional pole term 
which means that we also have to compute /_i(l). Therefore for z = 1 the phase 
space integrals ( |3.3CI| ) have to be computed even up to one order higher in powers 
of e than is needed for those which are integrable ai z = 1. Since fi^z) — fe{l) is 
integrable at 2; = 1 we can replace in ( |3.34| ) fe{z) by /(I) and one only has to consider 
the integral 

2m-l „i f (T\ 

T^-l = Y: / dzDf (^) (1 - .)f-^ (3.36) 
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which can be written as 



(m) 



2m-l 

E 



2 

+ —e 
m 



If we define the distribution (see [ ^ ) 
'ln^fl-2) 



V,{z) 



I- z 



by 



dz Vi{z)g{z) 



^0 I — z 



dz— '-{g{z)-g{l)), 



one can rewrite (|3.37|) in the following way 



(m) 



dz 



£=0 4=0^-' ■' 



(3.37) 



(3.38) 



(3.39) 



(3.40) 



where J-'xq stands for the soft gluon bremsstrahlung contribution which is given by 
(see the definition in 



^(r?i),soft 



5(1 



2 m- 1 



m 



-'Ml] 



(3.41) 



In order (m = 2) the highest order pole term which can occur in ( p.41| ) is repre- 
sented by 1/e^. The latter is cancelled by similar terms originating from the virtual 
gluon contributions given by the two-loop vertex corrections in fig. ^. Finally we 
want to emphasize that the type of singular terms in ( |3.34| ) only occur in JF^^ and 

T^^g and are absent in JF^^ (p = g, g) or J^k,g ^ {k = T, L). 

Adding all virtual-, soft- and hard-gluon contributions, the IR divergences cancel 
while computing the parton structure functions jF^p(z,Q^,£) which is in agreement 
with the Bloch-Nordsieck theorem. The left-over divergences are removed by coupling 
constant renormalization and the C-divergences are factorized out of J-'k^p{z,Q'^,e) 
leaving us with the coefficient functions which are finite in the limit e 0. These 
two procedures will be carried out in the next section. 



'''Notice that T^^^ — Ta.q = because of charge conjugation invariance of the strong interactions. 
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4 Determination of the coefficient functions in the 



MS- and the annihilation scheme (A-scheme) 



In this section we determine the coefficient functions of the process (|2.1| ) by applying 
coupling constant renormalization and mass factorization to the parton fragmenta- 
tion functions Tk,p (p = 9, g) which are computed up to order a\ in the last section. 
These coefficient functions have to satisfy renormalization group equations. One can 
formally solve these equations order by order in by writing the renormalization 
group functions like the beta-function (5{a^ and the anomalous dimension ^i^{a.^ 
{hi = 1)9) as a power series in ag. In this way one can algebraically express the 
coefficient functions into the coefficients of the power series. Using the mass factor- 
ization theorem which holds in every renormalizable field theory for all leading twist 
two contributions, one can also express the parton fragmentation functions Tk,p into 
the same coefficients. Our calculations described in the last section have to satisfy 
the algebraic expressions of j^k,p at least up to pole terms (l/e)"* which is a minimal 
requirement for the correctness of our results. 

Before presenting the algebraic expressions for J^^^p we have to decompose them ac- 
cording to the flavour symmetry group. Convoluting the parton structure tensor 
wjlu'^ ■* ( p.7| ) with the bare parton fragmentation densities D^{z) we obtain the fol- 
lowing functions 



P='l,q,9 9i>92=l 



<?1 "^92 



1 dz 



z P \z 



Dn-)PkA^,Q\e), {k = T,L), (4.1) 



P='l,Q 91:92=1 

^^AHm^^_^(^^g2^^)_ (4.2) 



The reason that we call Dp 'bare', originates from the fact that the C-divergence 
which are removed from JF^ p via mass factorization will be absorbed by Dp so that 
the latter are dressed up to the phenomenological fragmentation densities defined in 
( p.4|) , (|2.5|) . The hadronic fragmentation functions defined in ( p.21| ) are obtained by 
contracting the parton structure tensor W^iu''^ ^ (|3.7|), after convolution by Dp , with 
the leptonic tensor due to the subprocess e"^ + e~ V{V) where one also has to 
include the vector boson propagators given by Z(Q^)~^ in (|2.12|) . 
The contributions to J-'k,p {p = q, g) can be distinguished in a flavour singlet (S) and 
a flavour non-singlet (NS) part. Equations ( [4.1|) , ( [4.2|) can then be written as 



1 r "-f 



15 



qi = l p=l 131 = 1 p=l / 



n f . rif 

iji=i p=i 



+ (^) ) Q^ ^) + (f ) ^a^, Q\ e) 



ik = T, L), (4.3) 



^ dz 
z 



p=i 



+ (E E <^ + E «r E 

(ji=i p=i (}i=i p=i 



H 



Here we use the same notation as introduced above 
for p = u,d, . . .. Further we have the relations 



■p(r) 



{k = T,L-r = NS,PS), 



(4.4) 

|) where p = 1,2, ... ,nf stands 



(4.5) 



-f-NS 



-f-NS 



-f-PS 



A,q 



0, 



^,9 



0. 



(4.6) 



Relations ([4.5|) , (f4.6|) follow from charge conjugation invariance of the strong inter- 



actions. The parton fragmentation function JF^^ is called the purely singlet part for 
reasons we will explain below. 

The function JF^ ^ {k = T,L), describing process ( p.l| ) where the gluon is detected 
{p = 9)) receives contributions from the graphs in figs. H, H, ^ Since the gluon is a 
flavour singlet jF^ g belongs to the same representation. The quarks qi and q2 in (|4.1|), 
which are directly coupled to the vector bosons V and V' respectively, automatically 
belong to the inclusive state when p = g m reaction ( |3.1D so that the sums over qi, 
q2, and p in ([4.1| ) have to be separately performed. 

The non-singlet part JF^^ {k = T, L, A) describing process (|3.1|) where the quark or 
anti-quark is detected {p = q or p = q), is determined by the graphs in figs. except 
for the combinations C^, and AD, BC (see below). Notice that groups B and D 
only contribute when the anti-quarks qi and q2 are identical. In the case of the non- 
singlet contribution the quarks qi and q2 can be identified with p (i.e. p = qi = ^'2) 



so that the sums over qi, q2, and p in (|4.1|) are now connected. The above diagrams 



also contribute to in the case of k = T, L when they are projected on the singlet 
channel and the result is the same as the one obtained for JF^^ so that we can set 



k,, 



•r-NS 
k,q- 



The groups and in fig. R| only survive if they are projected on 



the singlet channel. This is because the detected quark p is only connected with the 
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vector bosons V and V via the exchange of a gluon which is a flavour singlet. To 
show this more explicitly we have drawn the cut graphs contributing to the parton 
structure tensor wjiu'^ ■* which originate from groups C and D in fig. ^. Because 
of the purely singlet nature the groups C and D only contribute to JF^^ and their 



p and qi, q2 which determines the factor of m 
non-singlet contribution which we will call ^^^^ ( 



contribution will be called ^k^- Like in the case of ^k,g the quarks qi, q2 belong 
to the inclusive state since p 7^ qi, p 7^ q2- Therefore one can separately sum over 

%.3\j . Finally we have a special 
{k = T, L) and fj). The 
latter originates from the combinations AD and BC in fig. ^ which only appear in the 
case when the anti-quarks pi and p2 are identical. The corresponding cut graphs are 
drawn in fig. |10|. If one removes the dashed line, which indicates the integration over 
the momenta cut by that line, one obtains a closed fermion loop. This fermion loop, 
to which are attached two gluons and one vector boson V {V'), has the same prop- 
erties as the triangular fermion loops inserted in the virtual diagrams of figs. ^, ^. 
In fig. |l^ we have taken the example that the vector boson V couples to the cut 
fermion loop via the quark qi whereas V couples to the detected quark p (see also 
( [4.3|) , ([4.4|) ). Like in the case of the triangle fermion loops in figs. only the axial 
vector current can couple to the cut fermion-loop which rules out y = 7 so that 

we have in (|4.3|) V = Z whereas in ( [4.4| ) 



,(7) 



only V = Z remains. Since a"p^ 
we either can get V = 01 V = Z. Only when the above condition is satisfied the 
parton fragmentation functions :P'j^^ {k = T, L) and P'2^g^ can contribute to Pj^^'^'^ 

^ and F^'^^ {V = 7, Z) respectively. If we now in addition sum in fig. |1^ over 

(see (OD) so 



all quark flavours qi belonging to one family one gets ^^^^^ ^ a^f ^ 

that in this case the above contributions due to ^'1^^, -^'aq" ^^^^ vanish. Since one 
has to sum over all members of one family anyhow in order to cancel the anomaly 
appearing in the triangle fermion-loops in figs, ^^[we will do the same for the graphs 
in fig. |TU[ Therefore we do not have to calculate j^jj^^ and ^^^^ and they will not be 
included in our phenomenological analysis in this paper. 

Summarizing the above the singlet fragmentation function JF^^ [k = T,L) receives 
two kinds of contributions and it can be written as 



k 



ik = T,L). 



(4.7) 



After having specified the various parts to the parton fragmentation functions we will 
now list them below. Starting with the non-singlet part the parton fragmentation 
function expanded in the bare coupling constant read as follows 



'NS,(2) 



L,g 



5 



q: 



2 r 



+ P(0) 

^ qg 



.(1) 



+ C 



NS,nid,(2) 
L,q 



+ C 



_NS,id,(2) 



L,q 



qq 



L,q 



(4.8) 



JF 



■NS,(2) 



T,q 



^] Si (91 



1 ( ip(0) 

£2 12 



p(0) 

qq 



qq 
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, p(l),NS^ -(1) p(0) ^ -(1) I , ^NS,nid,(2) NS,id,(2) „^ (1) 



The convolution symbol denoted by ® is defined by 

{f®9){z)=l dzi / rf2;2 5(z - ZiZ2)/(2i)5'(^2)- 
Jo Jo 



(4.9) 



(4.10) 



The second order DGLAP splitting functions denoted by P/^^^ (z, j = g, (?) are different 
for deep inelastic structure functions (spacelike process) and fragmentation functions 
(timelike process). For the latter case they have been calculated in ||55| , |58| . In 
order to solve the Altarelli-Parisi equations of the fragmentation densities it is 
convenient to split them into two parts which in the MS-scheme are given by 



160_ , , 16 176 161 + ^2 

Voiz) \ z In z 

9°^^ 9 9 3 1 -z 



4 32 
6{l-z){- + ^a2)) 



l + z' 
l-z 



In z (12 + 16 ln(l - z) - IQ In z) - 40(1 - z) 



- (28 + 12^) In z + 4(1 + z) In^ z + 5(1 - z) (3 - 24C(2) + 48C(3)) 



+ CaCf 



^36 919 

16C(2))Po(^) + 8(1 + z)C{2) + 8(1 + z)\nz + 



9 



9 



+ (4 . + f m.) + .(1 - .)(^ + f C(2) - 24C(3)) 



(4.11) 



P^t'^'\z) = {Cl-'-CAC,) 



l + z^ 
l + z 



; In^ z - 32 In z ln(l + z) - 32Li2(-z) 



- 16C(2) + 32(1 -z) + 16(1 + z)\yiz 



(4.12) 



where Li„(x) denote the polylogarithmic functions which can be found in |5^. The 
splitting function PjJ^'^^^ arises when the (anti) quarks pi and p2 in reaction 

( p.l| ) become identical and it is only determined by the interference terms AB and 
CD m fig. S. Like the splitting functions we have also decomposed the second order 

' ■ -NS,i ■ 



coefficients c^g'^^"* {k = T,L) into two parts i.e. c^^''^''^'^^'' and cll'^'"^'*'^''. The latter is 
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due to identical (anti) quark contributions and like Pqq it originates from combinations 



AB and CD in fig. 



All coefficients cj^^^ (i 



0, 1) are computed in the MS-scheme 



indicated by a bar and they show up in the perturbation series of the coefficient 



functions as we will see below. The coefficients a[,^^ are presented in ( |3.18|) , (|3rT^ ) and 
(3o is the lowest order coefficient in the beta-function defined by 



-2a., 



/3o 



ff 



4 



Ttn 



(4.13) 



3 3 

where now stands for the renormalized coupling (see below). The purely singlet 
contributions (see ( [4.7|) ) are given by 



PS, (2) 



L,g 



An 



1 J lp(o) 
^ 2 ''^ 



+ C 



^PS,(2) 



L,g 



(4.14) 



PS, (2) 



T,q 



i J ipps,(i) 
2 '^'^ 



+ lp(0) 

2 «^ 



,(1) 1 , .PS,(2) lp(0)^ (1) 



where a^P and presented in ( ^.201 ) and ( ^.21] ) respectively. The above ex 



.(1) 



(4.15) 



pressions are determined by the combinations (non-identical (anti-) quarks) or 
and (identical (anti-) quarks) in fig. |^. The timelike splitting function Pq,j^'^^^ can 
be inferred from M and it reads (MS-scheme) 



— - 128 + 6Az + —z^ + 16(1 + z) In^ z 
9z 9 ^ ' 



(80 



128 , 

+ lAAz + —z ) \nz 



(4.16) 



From ( [4.7| ) we can now also obtain the singlet parton fragmentation function 



S,(2) 



{k = T, L). Adding eqs. ( ^I8|) and ( [4. 14|) provides us with JF^'^ whereas the sum of 
eqs. ( ^^ ) and ( |4.15| ) leads to JF^'^ . In the same way we obtain from ( |4.11| ), ( [4.12| ) 



and ( [4.16| ) the singlet splitting function 



pS,(l) 
99 



pNS,(l) I pNS,(l) pPS,(l) 



Finally the order contributions to J^k,g become 



-r-{2) 



471 



S 



q: 



(1) 



(4.17) 



.(1) 
-L,g 



(4.18) 
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(2) 



if 



(4.19) 



where the timehke sphtting function P^^'* in the MS-scheme can be found in BS 
It is given by 



p{i) _ 



4 + 36z + (-64 + Az) \nz + IQz ln(l - z) + {% - Az) In^ z 



+ (— - 16 + 8z) ln2(l - 2) + (— - 64 + 32^) In z ln(l - z) 

z z 



128 128 
+ ( 128 + 64z)Li2(l -z) + ( + 128 - 64z)C(2) 



+ CaCf 



136 „ 352 2 , 48 64 2^ , 

h 40 - 8^ + ( h 64 + 72;z H z^) In ;z 

9;z <d ^ z 3 ^ 



- 16z ln(l - 2) - (— + 16 + 24z) \yi^ z+{ + 16 - 82) ln^(l - z) 

z z 

32 128 

+ ( + 32 - IQz) In z ln(l - ^) + ( + 128 - 64z)Li2(l - z) 

z z 

32 32 128 
+ (— + 32 + 16z)Li2(-;z) + (— + 32 + IQz) In z ln(l + z) + ( 



- 96 + 642)C(2) 



(4.20) 



The pole terms {l/e)'"^ showing up in the parton fragmentation functions !Fk,p {k = 
T,L] p = q, g) are due to UV and C-divergences. In order to get the coefficient func- 



tions corresponding to the fragmentation process ( |2.1| ) these singularities have to be 
removed via coupling constant renormalization and mass factorization. 
The coupling constant renormalization can be achieved by replacing the bare (un- 
renormalized) coupling constant dtg by 



47r 



A-n 



An e V A* / 



2\ 



(4.21) 



where R represents the renormalization scale. After having removed the UV singu- 
larities the remaining pole terms can be attributed to final state collinear divergence 
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only because jFk,p is a semi-inclusive quantity. The latter singularities are removed 
by mass factorization which proceeds in the following way 



k,q 



■J k,q 
^k,g 



Qq 



2T 



(pNS 
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gq 



gg 



C 



k,g, 



(4.22) 
(4.23) 
(4.24) 



with r 



gq 



r - r 

^ gqi ^ qg 



Tqg. The quantities Tij are called transition functions in which 
all C-divergences are absorbed so that the fragmentation coefficient function C^^p are 
finite. Both functions are expanded in the renormalized coupling constant as{R^) and 
depend explicitly on the renormalization scale R and the factorization scale M which 
implies that they are scheme dependent. If we expand Tij in the unrenormalized 
coupling constant dg the expressions become very simple. Choosing the MS-scheme 
they take the following form 



-NS 

qq 



ip(o) 



+ 



qq 



qq 



gq 



qq 



y^qq 



[Ml 

47r ' V 



M2 

e/2 



+ 



1p(o) 



ipNS,(l) 1 pNS,{l) 

2 11 2 



+ 



a 



£1 



M2 

7^ 



(4.25) 
(4.26) 
(4.27) 



99 



99 



1p(i) 

2 



(4.28) 



99 



47r^^ 





r 




V7^. 




_2£ 



(4.29) 



99 



£/2 



1p(o) 

^^99 



(4.30) 



where the 1 in ( ^.25| ) and ( [4.3(]| ) is a shorthand notation for 5(1 — z). Notice that 
we have expanded the Tij above in sufficiently higher order of in order to get the 
coefficient functions finite. Therefore the computation of J-'k,p allows us to determine 
the DGLAP-splitting functions P^g^'^^\ Pgg^'^^\ Pqq^''-^\ and Pgq' in an alternative way 
which is different from the method used in [B3|, BSI. Further the transition functions 
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satisfy the following relations which originate from energy momentum conservation 



dzz{Tl{z) + Tg^{z)) 



! dzz{Tgg{z)+2nfTqg{z)) = l. 
Jo 



(4.31) 



(4.32) 



If we substitute JF^ p ( [4.22| ) - ( [4.24[ ) into eq. ( [4. 3D the C-singularities are absorbed by 
the bare fragmentation densities as follows 



■pNS 



D 



H 



r 



qq 



gq 



D 



D 



H 



H 

9 ' 



(4.33) 
(4.34) 
(4.35) 



Here D^^p and denote the non-singlet and singlet combinations of the parton 
fragmentation densities as defined in (|2.7|) , ( p.6|) . The same definition holds for the 
bare densities -D^g^ and . The densities -D^Sp' -^s^' depend on the renor- 

malization scale R and the mass factorization scale M which are usually set to be 
equal. 

Substituting eqs. (|42^)-(|4:23|), (|4:33|)- (|4:35|) in and using (|3) we obtain after 



rearranging terms the structure function f'^'^ \x,Q'^) expressed into the renormal- 
ized parton fragmentation densities and the fragmentation coefficient functions 
{p = q,g). 



^ dz 

z 



>- p=i 



X 



1 \ z / 

p=i 



(4.36) 



where we have chosen R = M. 

Like the parton fragmentation functions J^k,p in eqs. ( [4.8D -( ^?TU ) we can express the 
coefficient functions C^^p {p = q, g) into the renormalization group coefficients. In 
the MS-scheme they take the following form. The non-singlet coefficient functions 
become 



-pNS _ as 
^^■^ - 4^ 



.(1) 



+ 



Air 



1 



(1) 



Lm + c 



NS,(2),md 
L,q 



+ C 



-NS,(2),id 



L,g 



(4.37) 
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^Ns a. 



1 



tts 



-NS,(2),nid -NS,(2),id 



T,9 



The singlet coefficient functions are given by 



^-ps fas\ 



AttJ 



^P(O) ^ Jl) 1 . , -PS,(2) 



(4.38) 

(4.39) 
(4.40) 



+ c 



-PS, (2) 



The gluon coefficient functions become 



C 



L,g 



as 
An 



.(1) 



+ 



47r 



-Po-cil + lP^^cil 



,9 



(4.41) 



(4.42) 



C- 



47r 



+ 



47r 



P^XP^? + Pi^) 



+ c 



(2) 
T,3 



Further we have defined 



(4.43) 



(4.44) 



In the case M ^ R the resulting coefficient functions can be very easily derived from 
the above expressions ([4.37|) - ([4.43|) by replacing 



(4.45) 
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The explicit expressions for the coefficient functions ( [4.37| )-( ^^ ) are hsted in ap- 
pendix A. 

Besides the MS-scheme one also can compute the coefficient functions in the so called 
annihilation scheme (A-scheme) It is defined in such a way that F^/i?ee (see 

( f^.20D , ( |2.22| )) does not get any corrections at = = Q'^. In the A-scheme 
the transition functions Tij are related to the ones in the MS-scheme denoted by Tij 
(see (|4:25| ) - (|4:30|) ) as follows 



99 



where Z^q, are given by (see eqs. (2.61), (2.62) in [|^) 



7NS 
^99 



p— 1 (P"^ p— 1 (P 
^ee ^q ^ee ^9 



(4.46) 

(4.47) 
(4.48) 







The coefficient functions (r = NS, S, i = q,g) correspond to the structure func- 
tion defined in ( p.22| ) and they are given by 



C 



(4.49) 



The coefficient functions in the A-scheme, denoted by Cfc_p, are related to the ones 
presented in the MS-scheme in the following way {k = T, L) 



^k,q — \^qq ) ^k,q^ 



'Ck,i — {Z ^).. CfcJ. 



(4.50) 
(4.51) 



Expanding all coefficient functions and Ree in ag the former take the following form 
in the A-scheme 



-L,q 



as 
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^'1 471- 



(4.52) 
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(4.53) 
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iftpi? } il, + { Pi? - A4'i + ip,?' ® 4'i + Pi? ® 41 



-ip,l,?'®ci')U„-cS + 4"®41, 



(4.57) 



where 1 is given by ([4.44|) and the coefficients R^^^ show up in the perturbation series 
for i W^ - 

«» = l + S«'"+(S)'«'^'- 

Notice that we have expressed the above coefficient functions into the renormahzation 
group coefficients Pij \ Ck \ presented in the MS-scheme. From ( [4. 5 2| )- ( ^7571 ) we infer 
that at = {Lm = 0) the coefficient functions in (|4.49| ) become 

Cf=li?ee, Cf = 0, C, = 0. (4.59) 
In any scheme the coefficient functions satisfy the renormahzation group equations 



d d 



DM 



da. 



■-,{m) 



(4.60) 



with k = T,L and i,j = q, g. Further we have defined the MeUin transforms 

d^\QyM^)= f dzz^-'Ck,i{z,QyM'), (4.61) 

- [ dzz"^-^ Pij{z), (4.62) 
Jo 
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and introduced the following the notations 



k,q 



■'k,g 



2 



(4.63) 



The quantities 7^™'' are the anomalous dimensions corresponding with the timelike 
cut vertex operators of spin m. Like the timelike splitting functions Pij they are 
scheme dependent. The relations between the anomalous dimensions obtained from 
different schemes can e.g. be found in eqs. (3.82)-(3.86) in [ p^ . 
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5 Results 



In this section we will discuss the order contributions to the longitudinal and 
transverse cross sections and their corresponding fragmentation functions. In partic- 
ular we investigate how the leading order (LO) longitudinal quantities, which already 
exist in the literature [|^, [Q, are modified by including the order contribu- 



tions. We will do the same for the transverse quantities for which a next-to-leading 
order (NLO) result already exists. Further we study the dependence of the above 
quantities on the mass factorization scale M and the renormalization scale R and 
show that the sensitivity to these scales becomes less when higher order corrections 
are included. Before we proceed we want to emphasize that with all higher order 
QCD corrections at hand it is only possible to perform a full NLO analysis on the 
cross sections and the fragmentation functions. The order contributions also allow 
for a next-to-next-to-leading order (NNLO) analysis of the transverse cross section 
c"r(Q^) but not for the transverse fragmentation function For the latter 

one also needs the three-loop timelike splitting functions which have not been calcu- 
lated yet. Therefore the order contributions to Ft{x,Q^) have to be considered 
as an estimate of the NNLO result. Nevertheless we will use the notation F,^^^^ to 
indicate the order a1 corrected transverse structure function. 

The longitudinal and transverse cross section a^iQ"^) {k = T, L) defined in ( ^.181 ) are 
obtained from the coefficient functions calculated in the previous sections as follows 



The results are 
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+ nfCFTf\ 16C(3) + 



(5.3) 



Addition of ctl and ar yields the well known answer crtot(<5^) (see ( p.l9| ) and (|2.20|) ) 
which is in agreement with the literature ||3^ (see also ^^). Hence ( W^ and i f>.3\j 
provides us with a check on our calculation of the longitudinal and transverse coef- 
ficient functions. Notice that in lowest order crtot(Q^) only receives a contribution 
from the transverse cross section ( [5.3| ) whereas the order contribution can be only 
attributed to the longitudinal part in (|5.2| ). In order both ctl and ax contribute 

to CTtof 

Because of the high sensitivity of expression ( |5.2|) to the value of as, the longitudi- 
nal cross section provides us with an excellent tool to measure the running coupling 
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constant. 

To illustrate the dependence of the cross sections on the running coupling constant 
we have plotted the ratios 

2 



Rl{Q' 



0"tot(Q^ 



Arc 
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+ CaCf 



lllng-44C(3) + ip 



+ UfCpTf { 4 In ^ - 
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74 

y 



and 



0"tot(<5^ 



<yL{Q\ 
<7tot(<5^ 



(5.4) 



(5.5) 



as a function of Q (CM-energy of the e^e~ system) in fig. |TT] and fig. 12 respectively. In 
fig. we have shown corrected up to order (Rjf) and Rl corrected up to order 
{R^^'~'). For R'l^ we have used as input the leading log running coupling constant 

a^^{nf,A^^l\R^) with nj = 5 and Ag^ = 108 MeV {a^^{Mz) = 0.122). The input 
parameters of R^^^ are given by the next-to-leading log running coupling constant 

5 and Ag^ = 227 MeV ia^^^iMz) = 0.118). Further 



we have shown the variation of i?/, on the renormalization scale R by choosing the 
values R = Q/2, Q, 2Q. Fig. |Tl] reveals that the order corrections are appreciable 
and they vary from 48% {Q = 10 GeV) down to 28% (Q = 200 GeV) with respect 
to the LO result. Furthermore one observes an improvement of the renormalization 
scale dependence while going from R^'^ to R^^^. In fig. |T2| we have plotted Rt ( P^D 
up to first order (i?^-^*^) and up to second order {R^^^"^) in the running coupling 
constant. As input we have used for R^^^ and R^'^^^ the coupling constants a 



LL 



and a^^^ respectively. The features of fig. 12 are the same as those observed in 
fig. |ri|. In particular R^^^^ becomes less dependent on the renormalization scale as 
is shown for R'^^^ . In figs. |ll], |l^ we have also presented the values Rl and Rt at 
Q = Mz measured by the OPAL-experiment |^ which are given by 



Rr 



0.057 ±0.005, 



Rt = 0.943 ± 0.005. 



(5.6) 
(5.7) 



One observes a considerable improvement in the ratios R^ {k = T, L) when the order 
contributions are included. However there is still a little discrepancy between 
taken at R = Q = Mz, and the data. This can either be removed 



NNLO 
T ■> 



R^^^ and R 

by choosing a larger Af^ or by including the masses of the heavy quarks c and b in 
the calculation of the coefficient functions. Also a contribution of higher twist effects 
can maybe not neglected (see p|). 



We now want to investigate the effect of the order contributions to the longitudinal 
and transverse fragmentation functions Fl{x, Q^) and Ft{x, Q^) as defined in ( |2.21| ). 
For our analysis we have chosen the fragmentation density sets in 



which 

will be called BKKl and BKK2 respectively. The input parameters for as and the 
QCD scale A are the same as given below (|5.5| ) except for BKKl where one has 
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those given above for R^^ and R^'"^ respectively. However for both F^^^ (order 



A;;^ = 190 MeV. The definitions for F^^ and F^^° are the same as 

MS 



a 



corrected) and F^'^^^ (order a1 corrected) we use af^^^(5, A^, i?^). In 



the fragmentation densities D!^[z, M^) have been determined for H 



+ K~ by fitting the total fragmentation function F{x,Q^) = Y2h F^{x,Q^) 
with H = 7r±, K^, P, P to the e+e" data in the range 5.2 < Q < 91.2 GeV. 
Here the proton and anti-proton contributions to the fragmentation functions have 
been estimated like 



with 



F,^+^(x,Q2) = (l + /(x))Ff +--(x,Q2 



/(x) = 0.16, in [63], 



(5.8) 



(5.9) 



/(a;) = 0.195 - 1.35(a;-0.35)^ in 



Further we introduce the notation 
Fk{x,Q^) = Y,F^{x,Q^), 

H 



(5.10) 



(5.11) 



where we sum over H = tt^ 



K+, K-, P, P. 



Notice that f{x) in ( p.lO|) becomes negative when x > 0.73 so that (x, 



ceases to be valid above this x-value. In ||6^ (BKKl) the fit has been only made to 
the TPC/27-data g (Q = 29 GeV) whereas in (BKK2) one also included the 
data coming from the ALEPH p8| and OPAL pO[ collaboration. Since the range 



I and OPAL 

of Q-values covered by the BKK2 parametrization is larger than the one given by 
BKKl the scale evolution of the fragmentation densities turns out to be better when 
the BKK2-set |^ is chosen. However this improvement goes at the expense of the 
description of the longitudinal fragmentation function Fi{x,Q'^) as we will show be- 
low. For each set there exists a leading log and a next-to-leading log parametrization 
of Dp{z,M^) {p = q,g). The latter is presented in the MS-scheme so that we have 
to choose the corresponding coefficient functions in appendix A. Further we set the 
factorization scale M equal to the renormalization scale R. 
In fig. H we have plotted Fl{x, Q^) in LO and NLO at M = Q 
the results with the ALEPH H and OPAL 



data. 



Mz and compared 
Here we have chosen the 



BKKl-set because the BKK2-set leads to a much worse result. The latter already 
happens in LO as was noticed in where one had to choose a very small factor- 



ization scale. Here we have chosen the BKKl-set. We observe that F^'^ is below 
the data in particular in the small x-region. The agreement with the data becomes 
better when the order corrections are included although at very small x F^^^ is 
still smaller than the values given by experiment. In the case of the BKK2-set (not 
shown in the figure) one gets a result which is far below the data. This was already 
noticed in fig. 5 of where one had to choose a very small factorization scale M 
(M = 20 GeV) to bring Fj^^ in agreement with experiment. In NLO the discrepancy 
between F^^'^, in the case of BKK2, and the data becomes even larger which is due to 
the kaon contribution. It turns out that the convolution of Dp^~^^ {z, M^) with the 
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order contribution from the coefficient functions given in ( p.4| ), leads to a negative 
pK++K xhis example illustrates the importance of the measurement of Fl{x,Q'^) 
and the higher order corrections for the determination of the fragmentation densities. 
We have also shown the results for F^^'^ and F^^^^ at M = Q = Mz in fig. |T4 



using the BKKl-set. Both fragmentation functions agree with the data except that 
pNNLO gg^g ^ little bit worse at very small x. Furthermore F^^^ and F^^^^ hardly 
differ from each other which means that the order a1 corrections are small. We do not 
expect that this will change when the three-loop splitting functions are included. One 
also notices that Fi constitutes the smallest part of the total fragmentation function 



F = Ft + Fl which can be inferred from figs. |13|, [1^. This in particular holds at 
large x where Ft » Ft. Hence a fit of the fragmentation densities to the data of Ft 
is not sufficient to give a precise prediction for Fl and one has to include the data of 
the latter to provide us with better fragmentation densities. This in particular holds 
for [z, M^) in the small z-region. The order contribution to Ft will certainly 
change the parametrization of the gluon fragmentation density given by ALEPH in 
2g and OPAL in 



To illustrate the effect of the order contributions to the coefficient functions cal- 
culated in this paper at various e'^e~ collider energies we have studied the i^-factors 



= 'n.^orZ^ (5-12) 



In fig. |T3|we have plotted (^) for H = 7r++7r- at Q = 5.2, 10, 29, 35, 55, 91.2 GeV 
choosing the BKK2-set since the latter shows a better scale evolution. From fig. [1^ 
one infers that the corrections are large at small x where they vary between 2 {Q = 
5.2 GeV) and 1.4 {Q = 91.1 GeV). The corrections become smaller when x increases. 
A similar plot is made for KJ[, '^'^ in fig. |l6l Here the order corrections are much 



smaller than in the longitudinal case except at large x where they are of the same 
size. Furthermore at low x the order corrections become negative (i^J '^^ < 1) 
which is already revealed by fig. for Q = Mz- Again the largest correction occurs 
at smallest Q. This can be mainly attributed to the running coupling constant which 
becomes large when Q gets small. In fig. we investigate the dependence of Kt 
( ^.13|) on the specific set of fragmentation densities used. The same was done for 



Kl ( 5.12|) in |3^ where we compared the BKKl-set with the one in which is 
presented in the A-scheme. It turned out that is very sensitive to the chosen 
parametrization. Choosing the set in [^] and M = Q = Mz it turns out that Kl 



is mildly dependent on x. Using the same input a similar observation can be made 
for Kt which shows a constant behaviour over the whole x-region (see fig. |l^). On 
the other hand (see fig. 2 in |]53|) and Kt steeply rise when x tends to one if 



the BKKl or BKK2 sets are chosen. Hence we conclude that the i^-factors heavily 
depend on the chosen parametrization for the fragmentation densities. 
In the next figures we study the factorization scale dependence of the fragmentation 
functions and show the decrease in sensitivity on the scale choice for M when higher 
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order corrections are included. 

In fig. |l8|we have plotted F^^^{x,Q'^) at three different scales M = Q/2, Q, 2Q 
where Q = Mz- Like in fig. O we have chosen the BKKl-set since in this case we 



get agreement with the data. From fig. 18 one infers that the scale variation of Fl is 



small and that all scales describe the data rather well. To show the improvement in 
the scale dependence more clearly it is convenient to plot the following quantity 

... , max {Fl{\Q),FaQ), F[{2Q)) - min (F^ig), F^jQ), Fa2Q)) 
^ ' average (F,^(iQ), Fj:{Q), Fl{2Q)) 

(5.14) 

for r = LO, NLO, NNLO and k = T, L. 

In fig. IS one can see that A^^^ < IS^ as long as x > 0.1 which implies that Ff^^^ is 
less sensitive to the scale M than F^^ . The same observation is made for Ft (fig. ^0]) 
and At (fig. P). In fig. |0| we have plotted F^^° at the same scales as above. At 



all three scales the data are described very well. In addition we also show F^^^^ for 
M = Q = Mz- At very small x we found that F^^^^ < F^^^ for all three scales 
and F^^^"^ is also below the data. In fig. pi] we see again an improvement while 



going from LO to NNLO except when x < 0.1. The reason that for x < 0.1 the scale 
dependence of F^^'-' {k = T, L) is much larger than the one found in F^^ can be 
found in where is stated that the scale evolution of the fragmentation densities 



is only reliable in the region 0.1 < x < 0.9. 

Finally we also study the scale dependence of the fragmentation functions at a lower 
energy. As an example we take the total fragmentation function F^ with H = tt^+tt^ 
( p. 22 ) and investigate its behaviour for different choices of the factorization scale M 



where again M = Q/2, Q, 2Q. Contrary to the previous plots we have chosen the 
BKK2-set which range of validity is bounded by Q < 100 GeV and 0.1 < x < 0.8. 
Further we take Q = 29 GeV and compare the theoretical result with the TPC/27- 
data 1^ . In fig. we show F^'^^^ at three different scales. The scale variation is 
small and only noticeable at large x. The data are very well described by i^^-^^o 
the three different scales except at very small x where the BKK2 parametrization is 
not reliable anymore. The same holds for F^'^'^^^ which hardly differs from F^''^^^ 
so that even at lower energies the order corrections are very small. To show the 
improvement of the scale dependence in a better way we have plotted A^'^*^, 
and gg_ 23. A comparison with fig. 21 shows that there is essentially 



no difference between the A^ (r = LO, NLO, NNLO, k = T, L) taken at low 
(Q = 29 GeV) and high energies (Q = Mz = 91.2 GeV). 

Summarizing our paper we have computed the order contributions to the longi- 
tudinal and transverse coefficient functions. The effect of these contributions to the 
longitudinal and transverse cross sections are large which allow us for a better de- 
termination of the strong coupling constant as- The corrections to the longitudinal 
fragmentation function F^ are appreciable too which has important consequences 
for the determination of the gluon fragmentation density F)^{z, M^). 
Furthermore one can now make a full NLO analysis of F^ . A NNLO description of 
the transverse structure function is still not possible because of the missing three-loop 
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DGLAP splitting functions. However the order contributions from the coefficient 
functions indicate that probably the NNLO corrections are very small. 
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Appendix A The coefficient functions in the MS- 

scheme 



In this appendix we will present the explicit expressions for the coefficient functions of 
the fragmentation process in ( |2.1| ) which are calculated in section 4 in the MS-scheme. 
In order to make the presentation self contained we also give the order contribu- 
tions c^^p ( |4.37| )-( 4.43| ) which have already been presented in the literature ||D, 
The coefficient functions C^^p {p = q, g) will be expanded in the renormalized cou- 
pling constant ag = as{M'^) where we have chosen the renormalization scale R to be 
equal to the factorization scale M. If one wants to chose R different from M, q;s(M^) 
has to be expressed into ^^(i?^) following the prescription in ( [4.45|) . In this paper 
we will only present the expressions for the transverse coefficient functions since the 
longitudinal ones are already shown in [Q. However we will make an exception for 



the order corrections to the non-singlet part C^^. For future purposes we want 
to split it into a part due to identical quark contributions (AB and CD in fig. |]) 



represented by c^g'^^"*'"^ and a remaining part given by c^^^^^''^^^ (see ( [4.37|) and the 
discussion below (|4.12|) ). The same will be done for the transverse coefficient c^^'*'^^ 



in ( |4.38| ). The expression for the non-singlet coefficient is very long and we will split 
it up to the various contributions. First we have the soft plus virtual gluon contribu- 
tions which are represented by the distributions 6{1 — z) and T'o(-2) ( |3.38D . They are 
indicated by C^J^I . The remaining part which is integrable at z = 1 will be called 

C^^g\^ where H refers to hard gluon contributions although C^^]^ also originates 
from subprocesses with (anti) quarks in the final state (see fig. ^ except for and 
D^). Following this prescription the non-singlet coefficient function is constituted by 
the following parts 
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24 

+ In 2 ln(l + ;z)) + (-4 + Az + IQz^ ^^)C(2) + (12 + 8^ - Sz^ 

5 

12 2 . X, . X .24 118 42 24 
+ —z"^) In^ z + 16(1 - z) ln(l - 2) + (-^ + — + —z^) Inz 

ODD 



24 46 46 24 , 

\ z H z^ 

hz h h 5 



(A8) 



Notice that the identical quark contribution in ( [A7| ) and ( [A^ ) carries the colour factor 
The purely singlet coefficient function becomes 



16 . . 16 2 I r2 



CpTf^ 8(l + z)ln^ + — + 4-4^-y^^ \ Lj^ 



3 32 

+ CpTf { 16(1 + 2)(Li2(l -z)+\nz ln(l - z) + z + {— + 



-fz^)Hl-z) 



,64 „ 32 2,, 16 184 136 

8 - 40z z^) Inz \ z 

'3z 3 32 3 3 



64 



-PS, (2) 



+ —Z >Lm + Cp^g 



(A9) 



PS,(2) _^ rjn 



1{1 + z)(65i,2(l -z)- 2Li3(l -z) + 2 ln(l - 2)112(1 - z) 



+ 6 In zLi2(l -z) - 2C(2) In 2 + In z ln^(l - z) + 3 W 2 ln(l - z) 

+ ^ In^ 2 + (^ - 8 - 402 - ^2')(Li2(l -z)+ ln2 ln(l - ^)) + (1^ + 4 
6 32 3 32 

16 Ox, 2. X .64 16 o,, 2 .32 

- 42 - — 2^) ln2(l - ^) + ( 14 - 142 + — 2^) In^' 2 - (— + 32 + 322 

+ |22) (Li2(-2) + In z ln(l + ;^)) + (-^ - 8 - 242 + |2^)C(2) + (-g 



184 136 64 , , ,32 400 208 256 , 

\ z H 2^) ln(l - 2) - ( \ \ 2 H 2^) In 2 

3 3 3 ^^ ^^32 3 3 9 ^ 
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160 236 140 1024 



21 z 3 ^ 3 27 



(AlO) 



The gluonic coefficient function is equal to 



Ct,p - Cf 



a, 
An 



--8 + 4^ ^ Lm + (-- 8 + 4^)(ln(l-^) + 21n^) h8 

z \ z z 



16 

; 16 + 8^) ln(l - z) + {8-4z)lnz + 8-2z 



7^2 



+ 



— - 48 + 24^) (Lis (1 - ^) + In ^ ln(l - z)) + (— - 32 + 16z) 
z z 



48 

•(ln2(l -z)- 2C(2)) + (24 - 12z) In^ z + ( + 64 - 12z) ln(l - z) 

z 



g 

+ (-16 + 20z) In z H 20 + 8^ 

z 



+ CaCf 



16 16 1 24 

— - 16 + 8^) ln(l - ^) - (— + 16 + 16;^) hiz-^ 
' z z oz 



+ 32 + 4^ + —z^ 



64 

— + 48^)Li2(l -z)- (64 + 16z) 
z 



32 96 
•lnzln(l - ^) + (— + 32 + 16z)(Li2(-^) + ln^ln(l + z)) + (— 

z z 

16 64 

- 64 + 48z)C(2) + ( 16 + 8z) ln\l -z)-{— + A8 + 56z) In^ z 

z z 



/ 344 32 2,, , , , 400 

+ (-^ + 96 - 82; + -z^) ln(l -z) + (- — 



32 

64 + 40^ + In ;2 

o 



188 32 100 32 ^ 



Lm + Cj,g 



(All) 



p;{2) _ f^2 



16 



+ 32 + 16z)(4,Si,2(--z) + 41n(l + z)Li2(-^) + 21n^ln2(l + 



32 160 

- In^ z ln(l + ^) + 2C(2) In(l + z)) + ( 192 + 32^)Li3(-z) + ( 

z z 

+ 240 - 120^)5i 2(1 - ;2) - — Li3(l - ^) + (— - 32 + 16;^) ln(l - z) ■ 
' z z 

32 

•Li2(l - ;s) + (48 - 24^) ln^Li2(l - z) + { + 64 - 32^) ln^Li2(-^) 

z 
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AO 40 20 , ,44 22 o 48 , _ , 

+ (tT - ^ + 1^ (1 - + - ^'2) In ^ + (- - 40 + 20^) • 
oz 6 6 6 6 z 

32 16 

• In ;2 111^(1 - ;2) + ( 8 + 4;^) In^ ;2ln(l - ;2) + ( 48 + 24;^) • 

z z 

112 96 

•C(2) ln(l -z) + (-16 + 82)C(2) In^ + ( 208 + 402)C(3) + ( 

z z 

96 64 

+ 80 - 28^)Li2(l -z) + ( + 80 - 12^) In z ln(l - z) + ( 

z 5z^ 

- 64 - 96^ + ^z^) {U2{-z) +lnz ln(l + ^)) + (-— + 56 - 14^) • 
5 z 

• ln2(l -z) + (-32 + 83z - -z^) In^ ^ + (H^ - 112 - 44^ + — /)C(2) 

S z 5 

/8 , ,144 418 262 16 , 316 604 

+ (- - 24z) ln(l - z) + ( + z z^) laz+ 

^z ^ ^ ^ ^ 5z 5 5 5 ^ 5z 5 

154 16 2 

H z z^ 

5 5 



+ CaCf 



16 

— + 32 + 16^)(-25i,2(-^) - 2 ln(l + z)U2{-z) 
z 



32 

+ In zU2{-z) - In z W(l + z)- C(2) ln(l + z)) + (— + 32 + 16z) • 

z 

•(Lis ("Y^) - Lis (y^) +lnzln(l-z)ln(l + z)+ln(l-z) • 

48 240 48 

■U2{-z)) + (— + 32 + l&z)U:i{-z) + ( 160^)51,(1 - z) + { — 

z z z 

48 192 

+ 80 + 40^)Li3(l - ^) - (— + 16 + 40^) ln(l - z)U2{l - z) - ( + 32 

z z 

... M T- /. N / 8 8 4 , , 3,^ , ,160 88 
+ 144^) In zU2{l - ^) + (3; - 3 + 3^) In (1 ~ ^) ~ ^"37 + y 

1 94. s 4s 

+ —z) In^^z- (- + 24 + 12z)lnzln^(l-,z) - (— + 64 + 48z) • 
?> z z 

■ In^ zln(l - + (— + 32)C(2) ln(l - z) ^ (— - 96 + 128z)C(2) \nz 

z z 

+ (- + 48 + 2Az) In^ z ln(l + ^) + (^ + 40 + 12z)C(3) + (-^ + 32 
z z oz 
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32 n, , , , 496 32 , , , 

+ 24z + —z^)U2{l -z) + (-— + 96 + 16z + —z^) Inz ln(l - z) 



,80 



32 



+ + 80 + 562 + —z'){U2{-z) + ln;2ln(l + z)) + 
+ \n\l -z) + i-^ -1Q + 2Z- ^z') In^ z + 

o oZ o 

32 2x^/oN /356 236 4 32 2^ ^ ^ 
- y^')C(2) + - — - 3^ - y ^ ) 1^(1 - ^) + ( 



172 



48 -8z 



136 






+ 40z 




496 


772 


17^ 


3 



172 256 4438 928 ^ 

H ;z H z^) lnz + — 36 - 106^ z^ 

3 9^ 27^ 27 



(A12) 



Finally we also present the non-singlet longitudinal coefficient function (see also eq. 
(17) in H) 



^Ns as 



+ 



AttJ 



Cl< (81n(l -2) -41n2 + 2 + 42)L 



M 



22 



+ CACFi - —Lm \+nfCFTj{ -Lm \+c 



-NS,{2),nid -NS,(2),id 



+ C 



L,5 



(A13) 



c 



NS,(2),nid 
L,q 



Cl I 16S'i,2(l -z)- 32Li3(-z) + 16 In zU2{-z) - 16C(2) ln(l - z) 



- 12Li2(l - z)+A\yiz ln(l - z) + A In^il - ;z) + (-10 + 8z + Az^ + -z^) 

5 

In^ z + (24 - 16z - 8^2 - — ^^)C(2) + (-^ + IQz - 8z^ - —z^) ■ 

dz z d 



•(Li2(-z) + In^ ln(l + z)) + (14 + iz) ln(l - z) 



24 78 32 

5z 5 5 



16 2x, 24 211 86 16 o 

H z^) \nz-\ \ z H z^ 

5 ' 5z 5 5 5 



+ CaCf \ - 8^1,2(1 - ^) + 16Li3(-^) -8 In zLi2(-z) + 8C(2)ln(l-^) 



+ (-^ - - + 8;^ + 4^2 ^ ^^3)(Li2(-^) + Inz ln(l + z)) 
bz'' z 5 

9 4 2 ^^ 46 , , ,12 22 4 

+ (2 - 42; - 22;^ - -z')i\n' z - 2C(2)) - - ln(l - ^) + (_ + _ + 



-^^) In 2; h 



12 2017 254 



5z 45 



z z 

15 5 



40 



+nfCFTf I -(ln(l -z)+lnz)- — + -zj 



(AM) 



NS,(2),id 



325i,2(l -z)- IGLisi-z) + 32 ln(l + z)U2{-z) 



24 

+ 16C(2) ln(l + z) + 16lnz In^fl + z)-8\n'^z ln(l + z) - 16C(3) + 

16 16 

16 - lQz + 8z'^ - —z^){U2i-z)+\nz\n{l + z)) + (4 + 8z - 4z^ 

8 2 XX . 24 112 8 16 2x, 24 64 104 

+ - 2C(2)) + (-5; H- - - 5-- + y^^) In. + - + - - 



16 2 



(A15) 



7PS 



The remaining coefficient functions Cj^^ and C^^^ can be found in eqs. (19) and (20) 



respectively of . 
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Appendix B The coefficient functions in the anni- 
hilation scheme (A-scheme) 



A glance at the coefficient functions in the A-scheme as presented in ( [4.52| )-( p:.57| ) 
reveals that they are all expressed into the MS-representation of the DGLAP splitting 



functions P/^^ and the longitudinal coefficients {p = q,g). The longitudinal 



(fc) 



coefficient functions in the A-scheme are already given in eqs. (21)-(23) of and we 
only add those formulae needed to construct the transverse coefficient functions Ct,p 
in the latter scheme. For C^i ^ ( |4.53|) we have to compute the following convolutions 



y2 



12Vq{z) + 95(1 - 2) - 16 ln(l - z) 



+ 81n;z- 10- 14z 



and 



(p(i)5(l-;.)-c«) 



cl 



(Bl) 



20Li2(l - z) - 41n2(l - z) + 41nzln(l - z) 



41n2 z - 16C(2) + (10 - Az) ln(l - z) + (4 - 8z) In^ + 18 + 6z 



(B2) 



Note that apart from R^^\ R^"^^ the non logarithmic {Lm) contributions to C^^ ( 4.52 ) 
and ([4.53|) have an opposite sign. The same observation holds for the coefficient 
functions C^^ (|4.54| ), ( [4.55|) and <Ck,g (lill), (lil^)- Fo^' ^k,g we need the convolutions 



192 \ /64 

— - 160 + 80^1 Li2(l-z)+ f — - 64 + 32Z 



\n^{l-z) + 



12^ 



96 + 48Z lnzln(l - z) + (32 - lQz)\n^ z 



64 \ 80 

+ ( -— + 64 - 32^ j C(2) + (32 - 82) ln(l - z) + (-64 + 32z) In ^ - — 



and 



+ 96 - 16;z 



(B3) 



16Li2(l - z) + 16 + 161nzln(l - z) 



/32 \ 32 

ln(l-;z)+ ( — + 16-16;zj lnz + — - 56 + 24Z 



(B4) 
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Figure captions 

Fig. 1 Kinematics of the process e^e~ H + "X" . 

Fig. 2 Born contribution given by the subprocess V "g" + q. 

Fig. 3 One-loop correction to the subprocess V — > "g" + g. Graphs with external 
self energies are omitted since they do not contribute in the case of massless quarks. 

Fig. 4 Graphs contributing to the subprocess V "g" +q + g and V —>■ "(7" + g + g. 

Fig. 5 Two-loop corrections to the subprocess V "g" + g. Graphs with external 
self energies are omitted since they do not contribute in the case of massless quarks. 

Fig. 6 One-loop corrections to the subprocesses V — > "g" +q+g and V —>■ "(7" -|-g-|-g. 
Graphs with external self energies are omitted since they do not contribute in the 
case of massless quarks and gluons. 

Fig. 7 Graphs contributing to the subprocesses V — > "g" + q + g + g and V — > 
+q + q + g. 

Fig. 8 Graphs contributing to the subprocess V "g" -|- g(l) + q + g(2). The cross 
( X ) indicates that the process is exclusive with respect to the quark denoted by "g" . 
If g(l) 7^ g(2) only combinations A and C have to be considered. When g(l) = g(2) 
combinations B and D have to be added to A and C. 

Fig. 9 Cut diagrams obtained from the groups C and D in fig. ^ contributing to the 
process V "g" + g(l) + q + g(2). 

Fig. 10 Cut diagrams resulting from the combinations AD and BC in fig. |^ con- 
tributing to V ^ "g" -|- g(l) + g + g(2) in the case that g(l) = g(2). 

Fig. 11 The ratio Rl = crr/o'tot- Dotted lines: Rjf; solid lines: R^^^ . Lower curve: 
R = 2Q; middle curve: R = Q; upper curve: R = Q/2. The data point at Q = Mz 
is from OPAL [01. 



Fig. 12 The ratio Rl = (yr/otot- Dotted lines: R^^°] R^^^°. Lower curve: 
R = 2Q; middle curve: R = Q; upper curve: R = Q/2. The data point at Q = Mz 
is from OPAL H. 

Fig. 13 The longitudinal fragmentation function Fl{x,Q'^) at M = Q = Mz- Dot- 
ted hue: Ff ^; solid line: The data are from ALEPH [g^ and OPAL The 
fragmentation density set is BKKl [Q. 

Fig. 14 The transverse fragmentation function FtIxjQ"^) at M = Q = Mz- Dotted 
line: F^^^; solid line: F^^^^- The data are from ALEPH l2l and OPAL \M. The 



fragmentation density set is BKKl [Q. 

Fig. 15 The ratio = p^'^^^/F^'^^ with H = tx+ + tx' at M = Q for different 
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values of Q. Upper dotted line: Q = 5.2 GeV; solid line Q = 10 GeV; long dashed 
line: Q = 29 GeV; short dashed hne: Q = 35 GeV; lower dotted line: Q = 55 GeV; 
dashed-dotted line: Q = 91.2 GeV. The fragmentation density set is BKK2 [63 . 



Fig. 16 The same as in fig. 19 but now for the ratio 
H = TT+ + n~. 



^H,NNLO I TpH,NLO 



with 



Fig. 17 The dependence of Kt on the fragmentation density sets dX M = Q = Mz 
Dotted curve: set from [4C]; solid curve: BKKl |63]; dashed curve: BKK2 ||64|| . 



Fig. 18 The mass factorization scale dependence of F^^^ aX Q = Mz- Lower dotted 
curve: M = 2Q; middle dotted curve: M = Q; upper dotted curve: M = Q/2. The 
data are from ALEPH [EH] and OPAL [R2|. The fragmentation density set is BKKl 
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Fig. 19 Sensitivity of F£ (r = LO, NLO) to the mass factorization scale represented 
by Al (|1D at Q = Mz- Dotted line: A^^; solid hne: A^^"^. The fragmentation 
density set is BKKl |3[. 



Fig. 20 The same as in fig. [IB] but now for F^^^ . Also presented is the curve for 
F^^^o (sohd hne). 



Fig. 21 The same as in fig. but now for A^- ( p.l4|) . Also shown is A;^^^*^ (dashed 
line) . 

Fig. 22 The mass factorization scale dependence of the total fragmentation function 
pH^NLO ^-^^Yi = 7r+ + TT" at Q = 29 GeV. Lower dotted curve: M = 2Q; middle 
dotted curve: M = Q; upper dotted curve: M = Q/2; sohd line: F^'^^^^ . The 
data are from TPC/27 [Q. The fragmentation density set is BKK2 [|6^ . 

Fig. 23 Sensitivity of F^'^ (r = LO, NLO, NNLO) to the mass factorization scale 
represented by A"^"" with = 7r+ + tt" at Q = 29 GeV. Dotted line: A^-^*^; solid 
line: A^'^^'^; dashed line: A"'^^^°. The fragmentation density set is BKK2 [g. 
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Table 1 

List of parton subprocesses in e~^e~ annihilation up to order a^. 



figure 



Parton subprocesses 











V ^ Q + Q 
a]: q + q 

V q + q+'Y 
al: V q + q 

V + 9 

q + q+'Y 

q + q+'Y+9 

'Y+q + q' + q' 



(one loop correction) 



(two loop correction) 
(one loop correction) 

(one loop correction) 
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Fig. 4 
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Fig. 5 
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Fig. 8 
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Fig. 9 
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Fig. 15 
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Fig. 16 
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Fig. 17 
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Fig. 20 
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